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Abstract 

In this paper we explicitly calculate fundamental solutions for the 
Tricomi operator, relative to an arbitrary point in the plane, and show 
that all such fundamental solutions originate from the hypergeometric 
function F(l/6, 1/6; 1;£) that is obtained when we look for homoge- 
neous solutions to the reduced hyperbolic Tricomi equation. 

1 Introduction 

The Tricomi operator 

d 2 d 2 



dx 2 dy 2 ' 

one of the simplest examples of a partial differential operator of mixed type, 
is: (i) elliptic in the upper half plane (y > 0); parabolic along the x-axis 
(y — 0); and (iii) hyperbolic in the lower half plane (y < 0). 

Our aim is to obtain explicit solutions in the sense of distributions or 
generalized functions of the equation 

(1.2) TE = 5(x-x ,y-y ), 

where S(x — x ,y — y ) is the Dirac function at (x , yo) an arbitrary point in 
the plane. A solution E of ( |1.2|) is said to be a fundamental solution relative 
to the point (xo,yo). 
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In a previous paper [fj we considered the case when x = y = and 
proved the existence of two remarkable fundamental solutions that clearly 
reflect the fact that the operator changes type across the x-axis. In this 
paper we study the general case when x = y ^ and compare our results 
to those of @. 

For sake of completeness and in order to make the reading of this paper 
independent of that of @, we briefly review the contents of that paper. It 
is known that the equation 9x 2 + 4y s = defines the two characteristics for 
the Tricomi operator that emanate from the origin. These characteristics, 
which are tangent to the y-axis at the origin, divide the plane in two disjoint 
regions D + and D_ (Figure 1) defined as follows: 

D+ = {(x, y) eM 2 : 9x 2 + % 3 > 0}, 
the region "outside" the characteristics, and 

D_ = {(x, y) e M 2 : 9x 2 + Ay 3 < 0}, 
the region "inside" the characteristics. 




Figure 1 

The first fundamental solution is defined by 

r c+(9x 2 + v)- 1 / 6 

(1.3) F+{x,y) = \ 

( o 

2 



in D+ 
elsewhere 



with 

< l4 > c +=-2wW44 ;i;i; 

and the second one is defined by 



;i.5) F_(x,y) 



C-\9x 2 + 4y 3 \- 1 / 6 in£L 
elsewhere 



with 

In the expressions of both C + and C_, the constant 

F(U;i ; i)- r(2/3) 



•6' 6' ' ' P(5/6) 

is the value of the hypergeometric function F(l/6, 1/6; 1; £) at £ = 1. 

We remark that in the constants C + and C_ were respectively denoted 

r(l/6) ^ 31X4/3) 



3-2 2 / 3 7r 1 /2r(2/3) * 2 2 /% 1 /2r(5/6)' 

At that time, we were unaware of the role played by the hypergeometric func- 
tion F(l/6, 1/6; 1; () in the study of fundamental solutions for the Tricomi 
operator. 

We also observe that according to formula ( |1.5|) , a perturbation at the 
origin spreads out to the entire region Z}_, inside the characteristics. On the 
other hand, according to formula (|1.3| ), the same perturbation spreads out 
to the whole elliptic region (y > 0) and also to the hyperbolic region (y < 0) 
outside the characteristics. Since the Tricomi operator is invariant under 
translations along the x-axis, the same phenomena take place for fundamental 
solutions that correspond to the Dirac measure 5 concentrated at an arbitrary 
point (a, 0) on the x-axis. 

It is well know that the Tricomi operator describes the transition from 
subsonic flow (elliptic region) to supersonic flow (hyperbolic region). In the 
extensive literature on the Tricomi operator one finds in the works of sev- 
eral authors, among others, Agmon [J] , Friedrichs , Gelfand [unpublished] , 
Germain and Bader M, Landau and Lifshitz B, Leray [O], and Morawetz 



12], a large number of examples of solutions to different problems that show 
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the interaction between these two regions. It seems that the two distributions 
F + (x,y) and F_(x,y) are the simplest of such examples. 

We now return to equation ( p..2|) . In view of the invariance of the Tricomi 
operator under translations parallel to the x-axis, the problem of solving that 
equation is equivalent to that of solving 

(1.7) TE = 6(x,y-b), 

where b is an arbitrary real number and d~(x,y — b) denotes the Dirac measure 
concentrated at the point (0,6). 

Consider the case b < 0, that is, the point (0, b) is situated in the hyper- 
bolic region. We introduce the characteristic coordinates 

(1.8) t = 3x + 2(-y) 3/2 and m = 3x - 2(-yf /2 , 

set £ = 2(-o) 3/2 , and to let a > be such that 3a = 2(-6) 3/2 . Note that 
(£o, —£q) represents the point (0, 6) in characteristic coordinates. As shown 




Figure 2 

in Figure 2, two characteristics with equations 

(1.9) 3(x-a) + 2(-y) 3/2 = and 3(x + a) - 2(-y) 3/2 = 

4 



pass through the point (0, b). In characteristic coordinates, these are the two 
half-lines 



(1.10) £ = £ , -oo < m < £ , and m = -£ , -£ < £ < +00 

that originate from the point (£q, —£q). 

The characteristics (|1.9| ) meet the respectively, at the points (a, 0) 

and (—a, 0). Two new characteristics originate from these two points, namely, 

(1.11) 3(x - a) - 2(-yf /2 = and 3(x + a) + 2(-yf /2 = 0, 

which we call reflected characteristics. The corresponding equations in char- 
acteristic coordinates are 

(1.12) m = £ , £ < £ < +00 and £ = -£ , -00 < m < -£ . 

The two characteristics through the point (0, b) and the reflected charac- 
teristics divide the plane into four disjoint regions denoted by Dj, D n , D IH , 
and Div, and illustrated in Figure 2. 

In this paper, we show existence of four fundamental solutions each sup- 
ported by the closure of the corresponding region. Of these solutions, only 
two have physical meaning: the one defined in the region Dj and the other 
defined in the region Dn. This situation is similar to what happens in the case 
of the wave operator in two dimensions where the two relevant fundamental 
solutions are the ones supported by the forward and backward light-cones. 
As we will see, it is the hypergeometric function F(l/6, 1/6; 1; () that plays 
a crucial role in defining these four fundamental solutions. 

The plan of this paper is as follows. We prove in Section ^ that the 
function 

(1.13) E(£,m-£ ,-£ ) = (£ + £^\£ Q -m)-^ F(\, 1 1; f ~ ^ + f ) 

6 6 (£ + £ ){m — £ ) 

is a solution of T h u = 0, where T h denotes the reduced hyperbolic Tricomi 
equation (|2.1|) . 

After replacing into ( |1.13| ), £ and m by their expressions in ( |1.8| ), we 



obtain the following function of x and y: 

(1.14) E(x,y;0,b) = e i7r/6 [9(x 2 - a 2 ) + Ay 3 - 12a(-y) 3/2 ]~ 1/6 x 

'1 1 9(x 2 -a 2 )+ V + 12a(-y) 3 / 2 \ 
X 6 ' 6 ' ' 9{x 2 - a 2 ) + % 3 - I2a(-yf/ 2 J 
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This is the function that generates the four fundamental solutions relative 
to the point (0,6). In Section |3] we define the distribution Ei(x,y;0,b) as 
the restriction of E(x, y; 0, 6), after multiplication by a suitable constant, to 
the region Dj and show that Ej is a fundamental solution of the Tricomi 
operator. The region Dj is entirely contained in the hyperbolic region where 



it is natural to use characteristic coordinates. Theorem |3.1| is then proved, 
via integration by parts and by using the results of Proposition [2.1| . As a 
consequence of it we show that, as (0,6) tends to (0,0), Ei(x,y; 0,6) tends, 
in the sense of distributions, to F_(x,y), the fundamental solution defined 

by (ED- 

The fundamental solutions supported by the closure of the regions D IH 
and Djv are defined and studied in Section [|. In both regions, it is neces- 
sary to take into account the singularities of E(x, y\ 0, 6) along the reflected 
characteristics. The results of Proposition [O , that describe the asymptotic 
behavior of the hypergeometric function F(l/6, 1/6; 1; () as |£| — > +oo, are 
then needed in the proof of Theorem [4.1| . We note that as 6 — ► 0, both 
fundamental solutions E in and E IV tend to 0. 

In Section |5], we define E n (x, y \ 0, 6) as the restriction of E(x, y; 0, 6), after 
multiplication by a suitable constant, to Djj and show that TEn = 5(0, b). 
Since En is complex valued and the Tricomi operator has real coefficients, 
both its complex conjugate and real part are also fundamental solutions 
relative to the point (0,6). The imaginary part of Ejj is then a solution of 
the homogeneous equation Tu = and we will call it a Tricomi harmonic 
function. 

Contrary to what was proved in Corollary |3.1j it is not true that the 
fundamental solution F + (x,y) defined by ( |1.3| ) is the limit of Ejj(x,y;0,b), 
or its real part, as (0,6) — > (0,0). It is necessary to take a suitable linear 
combination of En and its complex conjugate in order to achieve this result. 



As a final remark, we mention that in his paper [11], J. Leray described a 



general method, based upon the theory of analytic functions of several com- 
plex variables to find fundamental solutions for a class of hyperbolic linear 
differential operator with analytic coefficients. In particular, he showed how 
his method could be used to obtain, in the hyperbolic region, a fundamental 
solution for the Tricomi operator relative to a point (0, 6). He also produced 
an explicit formula for the fundamental solution in terms of the hypergeo- 
metric function F(l/6, 1/6; 1; £). Our method is simpler, more direct, and 
gives us global fundamental solutions that clearly reflect the change of type 
of the Tricomi operator across the x-axis. 
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2 A Special Solution to T^u = 

In characteristic coordinates (|1.8|) the Tricomi operator T becomes 

d 2 1/6 d d 



[2.1) T h 



dldm £ — m d£ dm ' 



and we call 7^ the reduced hyperbolic form of T. 

We now look for homogeneous solutions of the equation T^u = 0. Every 
homogeneous function of £ and m of degree A, a complex number, can be 
written as 

u(£,m) =l x (f)(t) 1 

where is a function of a single variable t = m/£. Direct substitution into 
(|2.1| ) shows that (p{t) must be a solution of the hypergeometric equation 

«(i - tw\t) + [(5 _ a) - (l - mm + ^(t) = o. 

As a solution of this equation, we choose the following hypergeometric func- 
tion F(—\, 1/6; 5/6 — A; t) extended, by analytic continuation, to the whole 
complex plane C minus the cut [1, oo). Since we are looking for fundamental 
solutions to the Tricomi operator we take for A the value —1/6, as previously 
indicated in our joint paper 0. Thus 

is a solution of T h u = 0. 

Let now {£q, mo) be an arbitrary point in M 2 and consider the change of 
variables 

i — mo m — mo 
" — * —„ t- rn — > — . 



t - £ m-V 
After unenlightening calculations, one can show that the function 

(2.2) E(i, m- to, m ) = (£-m )- 1 / 6 (4-m)-^F(i I l; ^"^?"7v 

6 6 (£ — m )(m — £ ) 
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is also a solution of equation (|2.1|) . This is the special solution to T h u = 
that we are looking for. 

Consider now the adjoint to equation (|2.1|) : 

_ d 2 v 1/6 dv dv 1/3 
^ ' hV ~ ~m^i + I^m K ~d£ ~!fr^~ (t-mf V ~ 

One can see that if v is a solution of equation fl2.3p , then u = (£ — m)~ l ^v 
is a solution of equation (|2.1|) . Since E(£, m; £ , m ) is a solution of ( |1.13| ), it 
follows that the function 

(2.4) R(£, m; £ , ™o) = (I ~ m) l/3 E(£, m; £ , m ) 



is a solution of the adjoint equation (|2.3|) . 

Proposition 2.1. R(£,m; £ ,m ) is the unique solution of T^v = that 

satisfies the following conditions: 
1/6 

(i) Ri = R along the line m = mo, 

£ — m 

-1/6 

(ii) R m = R along the line £ = £$; 

£ — m 

(Hi) R(£ ,m ;£o,rn ) = 1. 

Proof. Clearly, conditions (i), (ii), and (Hi) imply uniqueness for R. As 
m = mo, the argument of the hypergeometric function in ( |2.2| ) equals zero 
and so F(l/6, 1/6; 1; 0) = 1. Thus, along the line m = mo, we have that 



(2.5) R(£,m -£ ,m ) 



£-m \ 1/6 i fl W dt 
£ - m 



1/6 1/6 

where am = . Therefore Re = R along m = mo- In the same 

t — mo £ — m 

-1/6 

manner, one can see that R m = R along the line £ = £q. Finally, it is 

£ — m 

clear that R(£o, m ; £o, m ) = 1. □ 

Remark 1. The function R(£, m; £q, tuq) is the Riemann function of the 
operator Th relative to the point (£o,mo), (see 0). 



8 



If we go back to the notations in the Introduction, with (0, b) and b < 0, 
and let £ = 2(— 5) 3//2 , then m = — £q and we obtain from equation ( |2.2|) the 
equation ( |1.13| ) that, for further reference, we write as follows: 

(2.6) E{£, m; 4, -4) = + 4>)- 1/6 (*o - m)" 176 ^, ±; 1; C) 

o o 

where 

{27) c (£-£o)(m + £ ) 



(£ + £ )(m-£ )- 

Similarly, after replacing £ and m by their expressions in (|1.8j ) we obtain 



the function (|1.14 ) that we rewrite as 

(2.8) E(x, y; 0, 6) = e ^ e [9(x 2 -a 2 ) +% 3 - 12a(-y) 3 / 2 ]- 1/6 * F i; 1; c) 
where 

9(s 2 -a 2 )+4y 3 +12a(-y) 3 / 2 
{ ' } ^ 9(x 2 -a 2 )+% 3 -12a(-y) 3 / 2 ' 

We wish to analyze this function. Denote by r a the reflected characteristic 
3(x — a) — 2(— y) 3 / 2 = at (a, 0) and by r_ a the reflected characteristic 
3(x + a) + 2{~yf' 2 = at (-a, 0). 

Proposition 2.2. E(x, y; 0, b) is a real analytic function defined in R 2 \(r a U 
r_ a ). The singularities of E(x, y; 0, b) occur along the reflected characteristics 
r a and r_ a . 

Proof. 1. Let 

z = [3(x + a) + 2{-yf /2 } [3(x - a) - 2(-y) 3/2 ] = 9(x 2 - a 2 ) + % 3 - 12a(-y) 3/2 

be the denominator of ( in the expression Q2.9| ). For y < 0, z = on r a Ur_ a , 
it is negative in the region inside these characteristics, and it is positive in 
the region outside them. For y > 0, 2; is a complex number. 

It is possible to choose the argument of z so that it varies continuously 
in the region M 2 \(r a U r_ a ). Indeed, define (f) : K 2 \(r a U r_ a ) — > Si, the unit 
circumference, by 

_ 9(x 2 -a 2 ) +% 3 - I2a(-yf' 2 

< P\ x ^y) — j 
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where 

p = |g( x 2 _ a 2 ) + % 3 - \2a{-yf /2 \. 

Since IR 2 \(r a Ur_ a ) is contractible, <fi lifts to M, that is, there exists a contin- 
uous function 9(x,y) defined on IR 2 \(r a U r_ a ) so that 

9(x 2 - a 2 ) + % 3 - 12a(-y) 3/2 = pe"* 6 ^. 

We may take 8(x,y) to be equal zero outside the reflected characteristis and 
equal 7r inside them. It follows that 

2-1/6 = [9(x 2 _ q2) + Ay 3 _ l 2a (- 2/ )3/2]-l/6 

is well defined and real analytic in IR 2 \(r a U r_ a ). As for 

(\ 1 9(x 2 - a 2 ) + Ay 3 + 12a(-yf/ 2 \ 
V6'6' ' 9(x 2 - a 2 ) + % 3 - 12a(-y) 3 / 2 J 

it is also real analytic in the same region. □ 

The following remaks are in order: 

1) For y = we have 

E{x,0;0,b) = F{-,-,l;l)\x 2 ~a 2 \- 1/6 , Vx G (-a, a) 
6 6 

and 

E(x, 0; 0, b) = e^ 6 F(l ± 1; l)(x 2 - a 2 )- 1 ^ V x £ (-a, a). 

o o 

2) For ?/ > we have 

C=f = so |C| = 1. 

It is known (see f|) that the hypergeometric function F(l/6, 1/6; 1; £) is then 
given by its hypergeometric series and this series is absolutely convergent in 
the closed disk \(\ < 1. 

3) Along the characteristics 

3(x - a) + 2(-y) 3/2 = and 3(x + a) - 2(-y) 3/2 = 
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through the point (0, b) the function E(x, y; 0, b) is equal to 



E = 2~ 2 l\by)- l ' i 

and so it has a singularity of order — 1/4 at y = 0. 

4) Along both reflected characteristics r_ a and r a , the function E(x, y; 0, b) 
has a logarithmic singularity. This follows from Proposition |4.1|, in Section 
U that describes the asymptotic behavior of the hypergeometric function 
F(l/6,l/6;l;C)as|C|-+oo. 



3 The Fundamental Solution in the region Dj 



Let E(x, y; 0, b) be the function defined by the expressions (|2.8|) and ( p79| 
and define 



(3.1) Ej(x,y-0,b) 



1 : E(x, y; 0, b) inDj 



2 1 / 3 ' 

elsewhere. 



Since E(x, y; 0, b) is C°° in Dj and bounded on the boundary of Dj, it follows 
that Ej(x,y;0,b) is a locally integrable function and defines a distribution 
whose support is the closure of Dj. 

Theorem 3.1. Ej(x,y;0,b) is a fundamental solution for the Tricomi oper- 
ator T relative to the point (0, b). 

Proof. We must show that 

(E I ,T<p) = <p(0,b), G C C °°(M 2 _). 

Since E(x,y;0,b) is locally integrable in M. 2 , this is equivalent to showing 
that 



(3.2) // E I (x,y,;0,b)T(p(x,y)dxdy = <p(0,b),V<l>eC, 



oo /jn>2 



By introducing the characteristic coordinates i and m, noting that in the 
new variables the Tricomi operator ( [LID becomes 



ax^ a|/ z \_oiom I — m oz am 
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and that the Jacobian of the transformation is 

d(x,y) 1 



d(£,m) 2 1 / :i 3 2 (£-m) 1 / 3 ' 
we can see that the integral in ( |3.2|) is equal to 



(£ - m)V 3 E(e,rrv,e , -l Q ){^t m — r ! —<f>t+ T J —K) didm, 

i — m I — m 



where <p(£,m) is <p(x,y) in characteristic coordinates. We denote the last 
integral by / and write it as 



r-eo roo 1 /g 1 /g 

(3.3) / = -/ / R{e,m;e ,-e )(<p im - T ^<p t + T ^—<i> m )dedm, 

J-oo Jen Z-m £-m 



-oo J in 



m £ — m 

after recalling that (£ — m) 1//3 E(£, m; £ , —£ ) = R(£, m; £ , —£o 



-a,0) 




(a,0) 



,(0,h) 




Figure 3 

Next we perform several integrations by parts and take into account the 
properties of the function R(£, m; £q, m ), with m Q = —£ , as stated in Propo- 



sition 2A. First consider the term in Q3.3Q that contains cj)£ m . Integrating by 
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parts first with respect to m, say, and then with respect to £, one obtains 



(3.4) - / / R<j) em d£dm = (j)(£ ,-£o) 

+ - / — — (f){£, -£ )d£+- (j){£ ,m)dm 

6 J to £ + £q 6 £ -m 

-k) poo 

/ Rim4> d£ dm. 

oo J £o 

Second integrate by parts, relative to £, the term that contains <pi in 
and obtain 

i r—to r°° p 

(3.5) - / / (j) i d£dm = 



6 J-oo J to 2-m 
I ' ' ° -<fi(£o, m ) dm / [ ttj, \ t~ — )4>d£dm. 



6 J-oo £ -m 6 J-oo Ji d£\£-m / 

Finally integrate by parts, relative to m, the term that contains <p m in 
373) and get 



Y f—vo POO Jj> 

(3.6) --j ^ ^ —^dldm: 



6Je t + to ' 6 J-oo Ji dm\£-m / 

By adding ( |3.4j) , (|375|), and ( |3.6| ) we obtain that I = (p(£o,—£o) and this 
completes the proof. □ 

As a consequence of Theorem |3.1| we obtain the following result. 



Corollary 3.1. As (0,6) — > (0,0), the fundamental solution ( \3. 1\ ) converges 
in the sense of distributions to F-(x,y) the fundamental solution defined by 
fiTty and (\LQ). 

Proof. First note that whenever Re(c — a — b) > 0, the value F(a, b; c; 1) is 
finite and we have 

T{c — a)F(c — b) 
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Thus F(l/6,1/6;1;1) is well defined. As (0,6) -> (0,0), then %,-£ ) -> 
(0, 0), and so the limit of E(£, m; 4, -4) defined by (p|) and Q is 

On the other hand, at the limit the region Dj coincides with £)_, the region 
inside the two characteristics that originate from the origin, where Im = 
9x 2 + Ay 3 < 0. □ 

Remarks. 1) In our joint paper ||, the multiplicative constant appearing 
in the definition of F-.(x, y) was given by 

c 31X4/3) 



2 2 / 3 7r 1 /2r(5/6)' 



It is a matter of verification that this constant coincides with the one given 
by formula ( |1.6|) . 

2) Theorem |3.1| could have been proved in a different way, namely, by 



using the Green's formula ( |4.5| ) and by integrating along a suitable contour, 
as we did in our joint paper ||. This method which will be explained in the 
following section, is more adequate in proving existence of the fundamental 
solutions supported by the closure of the regions Du, Dm, and Diy. 



4 Fundamental solutions in the regions Dm 
and Diy 

For reasons of symmetry it is enough to consider only the region Dm. This is 
the region bounded by part of the characteristic 3(x + a) — 2(— y) 3 ^ 2 = 0, part 
of the characteristic 3(x — a) + 2(— y) 3 ^ 2 = and the reflected characteristic 
3(ac - a) - 2(-y) 3 / 2 = (See Figure 2). 

As we already pointed out, the function E(x, y; 0, b) is singular along the 
reflected characteristic 3(x — a) — 2(— y) 3 ^ 2 = 0. The nature of the singularity 
depends on the behavior of the hypergeometric function F(l/6, 1/6; 1; £), 
where ( = (£ — £ )(m + £ ) / (£ + £ )(m — £ ) , along the characteristic m = £q. 
This is revealed by the following result on the analytic continuation of the 
hypergeometric series F(a, a; c; Q (see Erdely p). The same result also gives 
us the asymptotic behavior of the analytic extension as \(\ — > oo. 
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Proposition 4.1. For ( e C with |arg(— £)| < n, we have 

(4.1) F(a, a; c; () = (-C)- B pog(-CMC) + w(0] 

where 

(AO) y(ft- r(C) V ^'^"^ -" 

{ } U[Q -T(a)T(c-a)^ o (n!)=» C ' 

wrat/i (a) n = T(a + n)/r(a), and 

(«) B (o= r . ™ -, ± iaU1 -; +a) " h nC ", 

T{a)T{c-a) ^ {n\) 2 

with h n = 2^f(l + n)-^f(a+n)-^f(c-a-n), #(() = r'(C)/r(C)- Moreover, 
both series u(() and v(() converge for |£| > 1. 

From this proposition it follows, as we will show in the proof of Theorem 
4~1~1 , that the singularity of E(x, y\ 0, b) along the reflected characteristics is 



logarithmic. We may now define 
(4.4) E HI (x,y;Q,b) = 



^E(x,y;0,b) in D m 



2V3- 

elsewhere, 



a distribution whose support is the closure of D in . We then have the follow- 
ing result 

Theorem 4.1. Em(x, y; 0, b) is a fundamental solution for the Tricomi op- 
erator T relative to (0,6). 

Proof. 1. We use the following Green's formula for the Tricomi operator 
(see§): 

(4.5) / / (ET(p-(pTE)dxdy = / E(y(p x dy-(p y dx)-(f(yE x dy- E y dx) 

J J D JC 

where D is a bounded domain with smooth boundary C. If TE = on D, 
then last formula becomes 



(4.6) 



/ / ETip dxdy = / E(yip x dy - ip y dx) - tp(yE x dy - E y dx) 

J JD JC 
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and the contour integral may be used to evaluate the double integral. Now 
throughout this section the domain D = Dm lies entirely in the hyperbolic 
region and so it is more convenient to express the countour integral on the 
right-hand side of ( |4.6|) in terms of characteristic coordinates. 
One can check that, in these coordinates, we have 



where 



yip x dy — i^ydx = A(i — m) 1 ' 3 ^ di — A(i — m) x l 3r i\) m dm, 



and A = l/2 2/3 . Similarly, 

yE x dy - E y dx = A(£ - m) 1/3 E £ di - A(£ - m) 1/3 E m dm. 
Thus the contour integral in ((Of), denoted by Ic, can be written as 



(4.7) I C = A {i-mfl\E^,-^E,)di-A \ {t-m) x l 3 {E^ m -^E m ) dm. 
Jc Jc 




> 




m = £q — s 



to + e 



2. We must show that 



Figure 4 



(4.8) (E, Tip) =lim // ET^dxdy = -2 1/3 <p(Q, b), G C C °°(R 2 ), 
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where D e is the intersection of an open disk that contains the support of tp 
and the region, defined in characteristic coordinates, as follows 

D n i,e = {{£, m) : £ > 4 + e, -£ + e < m < 4 ~ e}. 

Since (p = on the boundary of the disk, it follows from Green's formula, 
that we can replace the double integral in (|4.8|) by an integral of the form 
(|4.7| ), where the contour C consists of the oriented line segments A e , B t , and 
C e shown in Figure 4. 

We start by evaluating the simplest of these integrals, namely, the ones 
along B e and C e . 

3. The integral along B e . Denote by Ig e the integral along the line segment, 
£ = 4 + e with m varying from 4 — e to — 4 + e. Taking into account the 
orientation of B t , we obtain from formula (|4.7Q that 

/to— « 
(4 + e - m) 1/3 £(4 + e, m; 4, -4)V>m(4 + e, m) dm 
•■«0+e 

-A / (4 + e - m) 1/3 E m (£ + e, m; 4, + e, m) dm, 



where A = 1/2 2 / 3 . 

Integrating by parts the first integral in ( 4~9|) , we get 

(4.10) Jb £ = A(4 + e - m) l l 3 E{£, + e, to; 4, -4)V>(4 + e, m) 



3 



(4 + e - m)- 2/3 E(£ o + e, m; 4, —£o)ip{£o + e, to) dm 



/£ -e 
(4 + e - m) 1/3 E m (£ + e, to; 4, -4)^(4 + e, "i) dm 

From equation (|2.4|) it follows that 



1? 



i(£ _ m yy*E(£, to; 4, -4) + (I - m) 1/3 E m (£, to; 4, ~4)- 



3 

Substituting into (|4.10|) , we obtain 



(4.11) I S( = A(4 + e - m) 1 / 3 E(4 + e, to; 4, - W(4> + e, ™) 
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-2A 



[i2m(4 + e, m; « , -*o) + « 77-^ —\ J Wo + e, m) dm. 



6 (4 + e - m) 



Now, as e — > 0, the last integral tends to zero because, by Proposition 2.1 



Rm = —R/Q(£ — m) along the line £ = £q. On the other hand, the first term 
in ( 4. 11|) , namely, 



e -e 



A{£ + t-m) 1/3 E{£ + e, m; £q, —£o)4>{£o + e > m ) 

A(2e) 1/3 £(£ + e, £ - e; 4, + e, £ - e) - 

A(2£ ) 1/3 £(4 + e, -4 + e; 4, ~W(4> + e, ~4> + e) 

tends to -A(2£ ) 1/3 £(4, -4; 4, - W(4, -4), as e -> 0. Therefore, 



(4.12) 



lim Ii 



-1 

22/3 



1>(l 



0, -*oJ 



because £(£ , -£ ', £ , -£ ) = (2£ )~ 1/3 and A = 1/2 2 / 3 . 

4. The integral along C e . Denote by Iq e the integral along the line m 
—£o + e with £ + e < £ < +00. From formula (|4.7f ) we now obtain 

/■+00 

(4.13) I Ge =A (£ + £ - e) 1/3 E(£, -£ + e; £ , -£ )M^ -*o + e) d£ 

Jio+e 

r+00 

-A (£ + £ - e^Eeit, -£ + e; £ 0} ~W(£, -£ + e) d£, 

Jio+e 

As before, integrate by parts the first integral, use the formula 



R e = -(£- m)-V 3 E(£, m; £ , -£ ) + (£ - m) l/3 E e (£, m; £ , -£ ). 

to substitute for Eg in the second integral, and rewrite ( [4.13[ ) as 

-A{2£ Q y/*E{£ Q + e, -£ , -e; £ , -£ Q )^(£ + e, -£ + e) 



-2A 



£ +e 6 (€ + £ -e) 

As before, we obtain 



(4.14) 



lim/, 



Cc 



2 2/3 



-4) 
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5. The integral along A e . In this case, m = £q — e, £ varies from +00 to 
d + e, and the integral to be considered is 



-OG 



(4.15) Ia. = -A / {£-£ Q + e) l ' 3 E{£,£ -e-t Q ,-l )il) l (£,£ -e)d£ 

Jto+e 
r+00 

+A (£-£ + e) 1/3 E e (£, £ - e; £ , -£ )i/>(£, £ - e) all 

Our aim is to show that 

(4.16) ]imI Ae = 0. 

Once this is done, then by adding (|4.12| ), (|4.14j ), and (|4.16|) , we obtain ( ^8|) 
and the theorem will be proved. 

As we did in the previous two cases, we integrate by parts the first integral 
in (|4.15|) and get 

(4.17) I Ae = A{2t) l ' 3 E{£ + e,£ - e; £ Q , - W(4 + e,£ - e) 
A f +OQ 

+- (£-£0 + e)- 2 / 3 E(£, £ - e; £ , -£ )^(£, £ ~ e) d£ 
h +t 



3 

2A I (£-£ + e)^ 3 E t (l, £ - e; £ , ~W(£, £ - e) all 



In order to prove (|4.16|) we now must take into account the asymptotic 
behavior of both E(£, to; £q, —£q) and Ee(£, to; £0, —to) for m near £ and, more 
specifically, the behavior of the hypergeometric function F(l/6, 1/6; 1; £) and 
its derivative 

(4.18) F'(li;l;C) = ^F(^;2;C), 

as ( — ► 00. 

To see this, we start by expressing (|4.17 ) in terms of F(l/Q, 1/6; 1; () 



and its derivative. In order to simplify notations we write from now on, 
F(C) = F(1/6,1/6;1;C) and G(() = F(7/6, 7/6; 2; (). In the case we are 
considering m = £q — e so, in view of formula ( |2.7| ), we set 

_ (2l -e)(l-l ) _ <x(e) (2l - e)(l - l ) 

C " + " ~ Wlth " (£ + £ ) ' 
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Also, recalling formula ( j2.6|) , we have 

(4.19) E(£,£ -e;£o,-£o) = e~ 1/6 (£ + £or 1/6 F(^-), 

—e 

and so 



(4.20) 



,-1/6 



-(£ + £o)- 7/6 F( 



E e (£,£ -e;£ ,-£ ) = 

aje) e-V%(2£ -e) ^ 
— e 18 — e 



in view of ( 4. 18| ) and taking into account that 

d( = 2£ (2£ -e) 
d£ -e(£ + £ ) 2 ' 



By substituting ( |4.19|) and (f4.20|) into (^.17|) and combining the resulting 
integrals we obtain 



(4.21) I Ae = A(2e) 1/3 (2£ + e)- 1/6 e - 1/6 F( 



e_-2£o 
2£ + e 



)iP(£ + e,£ -e) 



+ Vo _ e)5 /a p ^(^/o + f 27 ; -i^Fj^m, ^ - e) di 



(£-£ o y/s(£ + £ o y e 



—e 



A £ (2£ e)'^ f + °° (l- £ ° + £ ) 1/3 ( a ^)^ C ( a{e) )H£ £ e)d£ 



As e — > 0, the first term on the right-hand side of ( [4.21D clearly tends to 
zero. We must prove that the difference of the two integrals in ( [4.21[ ) also 
tends to zero as e — > 0. For this, according to formulas (|4.1|) , ( [12"D , and ( |4.3|) , 
we have 



)1/6 F( 



cr e 



, /O-lej. /(Tie), .cr e 



— e 



— e 



and 



write them as a sum + jjj - log e • 1^ , where 



By using these two expressions we combine the two integrals in ( |4.21| ) and 
write t 

(4.22) 
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(4.23) if = 



<„+« C - («) r/e (( + «o) V "< 



and 



(4.24) if } 



^ , na .-i/6 ^_^±e^_ Me) 

„ ll+e (i-e Q y/*(e + e ) u (— )me °- e)di - 



Lemma 4.1. i?ac/i o/ £/ie integrals Pi , 1 < « < 3, is 0(e 1 ^ 6 ). 



r(i) 

converge for large values of \(\, we may write that 



Proof of lemma. Consider first the integral 1^ . Since the series in ( f4.2| 



= Wo + and U(^) = U + —-U(^), 

— e cr(e) — e — e cr(e) — e 

where w and U are bounded functions for small values of e. 

By substituting these two expressions into (|4.22|) , rewrite it as a sum of 
three terms: 

(4.25) if = 

A r + °° [f — fnA- f") -2 / 3 

{ 3 (24 - f)5/6 "» L (I-OTW,) '* 11 * 4 -^ 

^ /.+(» (p_p \ f \yz 
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N5/6 



?o + e) 



-2/3 



i 

9 



o(24- 



-1/6 



log(a(e))- 



o- e 



log(a(e))- 



%+e (£-£ y/ 6 (£ + £o 

The first term (inside the brackets) in ([4.25 ) tends to as e — > 0. Indeed, its 
limit as e — > is 



r A(2£ ) 5/6 M ^o(2£o)- 1/6 t/o 1 r (* - 4)" 5/6 , ,24(£ - £ ) 



+ £0 



log( 



+ 4 



3 9 
The integral on the right-hand side is finite. On the other hand 
1 . „ 1 6 



u 



r(i/6)r(5/6) 



and Un 



r(5/6)r(7/6) r(i/6)r(5/6)' 



are the constant terms of the series u(Q and U(Q given by ( |4.2|) . It is a 
matter of verification that the quantity inside the brackets is equal zero. 

Next, consider the second term in (|4.25|) . After the change of variables 
£ — £q = te, one can see that the absolute value of that term is bounded above 



by 



Ce i/6 



/+00 
(t + i)- 2/3 r 7/6 logtdt, 



with C a constant independent of e. With a similar calculation, one can show 
that the third term in fl4.25|) is bounded above by 



Ce l/6 



/+00 
(t + l)V3 r 13/6 \ Qgtdt} 



with C another constant. Therefore, 1^ 



The expression ( |4.24| ) for similar to that of I^, is even simpler 
because of the absence of the factor log(cx(e)). Thus, with a similar proof we 

Oie 1 / 6 ). Finally, the integral if] given by flOp is 

analogous to Ijp with u and U replaced by v and V the power series defined 
by ( [4.3| ). The proof of the lemma is then complete. □ 

implies 

□ 



(3) 

also conclude that I A 

r(3) 



Recalling the expression ( 4.22 ) for one sees that Lemma 
( |4. 16| ) which completes the proof of the theorem. 
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5 Fundamental Solutions in the region Djj 

In the same manner as we did in the previous sections, define the distribution 



(5.1) 



E n (x,y;0,b) 



Y^E(x, y; 0, b) in D H 



2V3 







elsewhere, 



whose support is the closure of the region Djj. Our aim is to prove the 
following result. 

Theorem 5.1. Eu(x,y;0,b) is a fundamental solution for the Tricomi op- 
erator relative to the point (0,6). 

Proof. The proof is, with few modifications, analogous to that of Theorem 



O. As before, it suffices to show that 



(5.2) (E, Tip) = lim / / ETp dx dy = 2 1/3 cp(0,b), Vy? G C C °°(I 



where the domain of integration D e is defined as follows. Let D be an open 



'■A 




Figure 5 
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disk centered at the origin and with radius sufficiently large so that it contains 
the points (—a, 0), (a, 0), and (0, 6) and the the support of <p. Let D H e be the 
region obtained from Du by shifting it by e along the characteristics and by 
removing two small half disks centered at the points (—a, 0), (a, 0), as shown 
in the Figure 5. D e is then the intersection of D and Du e . 

By virtue of the Green's formula (|4.6|) , the double integral in the expres- 
sion (|5.2|) is to be replaced by a contour integral along the oriented paths A' t , 
7g, B' e , B t , %, and A e . 

The integration along A e is similar to the one along A t calculated in 
Section [| and one can see that its limit, as e — > 0, is zero. The same is true 
for the integral along A' e . 

The integral along B e , similar to the integral along B 6 (Section |j), tends 
to ip(£ , -£ )/2 2/3 , as e -> 0. The same thing happens with the integral along 
B' e . The sum of these two limits is then 2 1//3 ^(£o, — £q) which is equal to the 
right-hand side of Recall that if>(£, m) = <p((£ + m)/6, -((£-m)/4) 2 / 3 ). 

To complete the proof one has to show that the limits, as e — > 0, of 
the integrals along 7^ and 7 £ are both zero. Since these two contours lie in 
the elliptic region of the Tricomi operator, it is more convenient to use the 
reduced elliptic form of the Tricomi operator, namely, 

T-— — — — 

e dx 2 ds 2 3s ds 

which one obtains from equation ( |1 . 1| ) via the change of variables x = x 
and s = 2y 3 / 2 /3, and the corresponding Green's formula in the variables x 
and s. It is then a matter of verification, which we leave to the reader, that 
along both contours the integrands remain bounded and consequently both 
integrals along 7^ and 7 e tend to zero. The proof is then complete. □ 

We observe that by virtue of Proposition [2.2| , the fundamental solution 
Eu(x,y;0,b) is a complex valued in the upper half plane (y > 0) plus the 
region in the lower half plane outside the reflected characteristics and it is 
real valued in the triangle with vertices (— a, 0), (a, 0), and (0,6). Since the 
Tricomi operator has real coefficients, the complex conjugate Ejj(x,y;0,b) 
as well as the real part of Ejj(x, y;0,b) are fundamental solutions for T. 

Contrary to what happened in Corollary [TT] of Section [| where we showed 
how the limit of Ej(x,y;0,b), as (0,6) — ► (0,0), tends to the fundamental 
solution F_(x, y) given by formulas (|1.5| ) and (|1.6| ), neither of the fundamen- 
tal solutions just obtained tend to the fundamental solution F + (x,y) given 
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by formulas ( |1.3| ) and ( |1.4| ). For this to happen we need a particular linear 
combination of En and En. Let A and fi be such that 

f Xef+fie-f =-l/3 1 / 2 
\ A +ti =1. 

Then, the following result holds: 

Corollary 5.1. Let E n = XEu + fiEn- Then E n is a fundamental solution 
relative to (0,6) that converges in the sense of distributions to F + (x,y), as 
(0,6) -(0,0). 

References 

[1] S. Agmon, Boundary Value Problems for Equations of Mixed Type, 
Atti del Convegno Internazionale sulle Equazioni alle Derivate Parziali, 
54-68, Triestre (1954). 

[2] J. Barros-Neto and I. M. Gelfand, Fundamental Solutions for the Tricomi 
Operator, Duke Math. Journal, 98, no. 3, (1999), 465-483. 

[3] R. Courant and D. Hilbert, Methods of Mathematical Physics, vol. I, 
Interscience Publishers, New York, 1953. 

[4] Darboux, Legons sur la Theorie des Surfaces, vol. II, Gauthier-Villars, 
Paris, 1914. 

[5] A. Erdely, Higher Transcendental Functions, The Bateman Manuscript 
Project, Vol. 1, 2, 3, McGraw-Hill, 1953. 

[6] K. O. Friedrichs, Symmetric Positive Linear Differential Equations, 
Comm. Pure Appl. Math., 11, 333-418 (1958). 

[7] I. M. Gelfand and G. E. Shilov, Generalized Functions, vol. 1, Academic 
Press, New York, 1964. 

[8] P. Germain and R. Bader, Solutions Elementaires de Certaines 
Equations aux Derivees Partielles du Type Mixte, Bull. Soc. Math, de 
France, 81, 1953, 145-174 (1953). 



25 



[9] L. D. Landau and E .M. Lifshitz, Fluid Mechanics; Course of Theoretical 
Physics, Vol. IV, Pergamon Press, London, 1959. 

[10] J. Leray, La Solution Unitaire d'un Operateur Differentiel Lineaire 
(Probleme de Cauchy, II), Bull. Soc. Math. France, 86, fasc. I, 75-96 
(1958). 

[11] J. Leray, Un Prolongement de la Transformation de Laplace qui Trans- 
forme la Solution Unitaire d'un Operateur Hyperbolique en sa Solution 
Elementaire (Probleme de Cauchy, IV), Bull. Soc. Math. France, 90, 
fasc. I, 39-156 (1962). 

[12] C. S. Morawetz, A weak solution for a system of equations of elliptic- 
hyperbolic type, Comm. Pure Appl. Math. 11 (1958), 315-331. 

[13] L. Schwartz, Theorie des Distributions, Hermann, Paris, 1950-51. 

[14] F. Tricomi, Sulle equazioni lineari alle derivate parziali di secondo or- 
dine, di tipo misto, Rendiconti, Atti dell'Accademia Nazionali dei Lincei, 
Serie 5, 14, (1923), 134-247. 



26 



